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p^- Abstract 



We study the optimization of a neutrino factory with respect to non-standard neutral current 
neutrino interactions, and compare the results to those obtained without non-standard 
interactions. We discuss the muon energy, baselines, and oscillation channels as degrees of 
freedom. Our conclusions are based on both analytical calculations and on a full numerical 
simulation of the neutrino factory setup proposed by the international design study (IDS- 
CN ■ NF). We consider all possible non-standard parameters, and include their complex phases. 



We identify the impact of the different parameters on the golden, silver, and disappearance 
channels. We come to the conclusion that, even in the presence of non-standard interactions, 
56 ! the performance of the neutrino factory hardly profits from a silver channel detector, unless 
^ I the muon energy is significantly increased compared to the IDS-NF setup. Apart from the 
dispensable silver channel detector, we demonstrate that the IDS-NF setup is close to optimal 
S^ . even if non-standard interactions are considered. We find that one very long baseline is a 
^ ! key component in the search for non-standard interactions, in particular for \e"^^\ and |e™ |. 
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1 Introduction 

In neutrino physics, three-flavor oscillations have been successfully used as a model explain- 
ing all relevant neutrino data, see, e.g., Ref. [1]. In particular, the solar and atmospheric 
oscillation parameters have been measured with very high precisions, and the reactor mix- 
ing angle ^13 has been strongly constrained. Future experiments will test this small angle 
further, and be sensitive to leptonic CP violation and the neutrino mass hierarchy (see 
Ref. [2] and references therein). The ultimate high precision instrument for these purposes 
might be a neutrino factory [3-5]. Using different baselines and oscillation channels, it can 
basically disentangle all of the remaining oscillation parameters [6-9] in spite of the presence 
of intrinsic correlations and degeneracies [5, 10-12]. Because of its high precision, it might 
be natural to ask how sensitive it is to non-standard physics. 

Non-standard neutrino interactions (NSI) could be such messengers of new physics beyond 
the Standard Model in the neutrino sector. In this work, we will focus on interactions of 
the form 

{I7^7''(l-7')^a}{/7p(l±7')/}, 

which may affect the neutrino propagation in matter [13-16]. Here, / is an electron or a first- 
generation quark. Such dimension six operators can be considered as an effective low-energy 
fingerprint of new physics at a higher energy scale, once the high energy degrees of freedom 
have been integrated out. Note that these operators involve the same in- and out-state 
charged fermion, which means that they produce neutral current-like interactions. Similarly, 
one can write down operators involving different in- and out-state charged fermions, leading 
to non-standard effects in the neutrino production and detection [17]. From the theory point 
of view, dimension six operators suffer from the problem that the neutrinos come together 
with their SU(2) counterparts in the Standard Model, which means that charged lepton 
flavor-violating (LFV) processes are introduced at tree level (unless the SU(2) breaking 
effects are large jll; see, e.g., Refs. [18,20-22]. For example, if a = e, /3 = r, and / = e, 
T decays into three electrons are a consequence, which can be strongly constrained by B- 
factories. However, this SU(2) relation can be avoided when dimension eight operators are 
taken into account above the electroweak scale [22,23]. In such a case, charged LFV effects 
appear only at one loop level, leading to less stringent bounds. Of course, there are ways to 
circumvent this reasoning, but so far there is no motivation to assume that the NSI should 
be large. Therefore, we focus on NSI constraints in this study. For a summary of current 
bounds from non-oscillation experiments, see Refs. [2,23,24]. 

NSI are also constrained from the current oscillation experiments. The effect of NSI in 
solar neutrino experiments has been studied in Refs. [25-27], and in atmospheric neutrino 
oscillation in Refs. [18,28-32]. Some numbers on the constraints can be found in these 
references. In addition, NSI have been discussed in the context of future neutrino oscillation 
experiments in Refs. [33-54]. They can be also tested in astrophysical neutrino sources, 
such as supernovae [55,56], and in the early universe [57]. The sensitivity of non-oscillation 



-'^ There is a possibility to construct the dimension six operator without SU(2) counter processes, assuming 
a charged SU(2) singlet mediation [18]. However, it is constrained by the measurement of the Fermi constant 
and the lepton universality [18, 19]. 



experiments to NSI, such as collider and neutrino scattering experiments, has been pointed 
out inRefs. [22,58-62]. 

As far as the neutrino factory is concerned, the problem of parameter correlations and 
degeneracies has to be taken into account if one parameter needs to be extracted from the 
information encoded in the event rates. NSI might therefore be confused with the standard 
oscillation parameters, and the sensitivity to the standard oscillation parameters might be 
affected in the presence of NSI [34,41]. The discovery reach to NSI (including production 
and detection effects) has been studied in Ref. [50] for a setup with only one detector. A 
major step forward was taken in Ref. [53], where it was pointed out that a detector at the 
magic baseline has an excellent sensitivity to interactions of the form t'e + / ^^ '^r + /• On 
the other hand, the silver channel v^ —>■ Vt [6,7] is very sensitive to this particular interaction 
at high energies [46] , because it is the leading order effect in that channel [2] . Therefore, it 
is yet unclear what the best strategy to measure specific NSI actually is, and what the role 
of the silver channel for NSI could be. We will clarify the contribution of different channels 
to specific NSI in neutrino propagation in this study both analytically and numerically. In 
addition, we discuss the muon energy and baseline optimization of a neutrino factory for 
NSI, and compare it to the optimization for the standard oscillation (SO) parameters. We 
fully take into account complex phases, and we include all possible NSI in the discussion - 
we will comment on specific parameters in the following section. We adopt the point of 
view that NSI must be small, since there is not yet any theoretical motivation for large 
NSI. Therefore, we only discuss sensitivity limits, and no discovery reaches. We focus on 
neutrino propagation effects for the sake of simplicity, which means that we assume that 
production and detection effects are either not present, or constrained otherwise (such as by 
a near detector). Our starting point for the optimization will be the current baseline setup 
for the international design study of a neutrino factory (IDS-NF) [2,63], which is designed 
for optimal discovery reaches for sin^2^i3, the neutrino mass hierarchy, and leptonic CP 
violation. 

Our study is organized as follows: We discuss all possible NSI in neutrino propagation in 
Sec. |2l identify the ones relevant for a neutrino factory, and describe the channels provid- 
ing the main sensitivities analytically. In Sec. [3l we introduce our performance indicators 
and simulation details. The impact of different channels on specific types of non-standard 
parameters is, in a full simulation, illustrated in Sec. HI We then study the impact of the 
muon energy in Sec. 0, where we put special emphasis on the silver channel. The base- 
line optimization of the two-baseline setup is, for standard oscillation physics, revisited in 
Sec. El and it is compared to the same optimization for the NSI parameters. In Sec. [TJ we 
summarize the sensitivities expected from a neutrino factory, and we conclude. 

2 Phenomenology 

We focus on non-standard propagation effects in neutrino oscillations. These can be phe- 
nomenologically described by neutral current-type non-standard interactions (NSI) 

£nsi = {GF/V2){eg){VpYLvo.}{h,Pf] 



+ {Gp/V2){effj*{V^YLup}{h,Pf], (1) 

affecting the neutrino propagation in matter. Here, e^^ = (e^^f)*, P G {L,_R}, L = 1 — 7^, 
R = 1 + 7^, and / stands for all possible fermions in Earth matter {u quarks, d quarks, 
electrons). This definition includes the possibility of different non-standard interactions 
with quarks and leptons, and different interactions for left- and right-handed couplings to 
the fermionso Note that, in general, e^^ are complex numbers for a 7^ /?, and real numbers 
for a = [5. Since there are about two nucleons (a proton and a neutron) per electron in 
Earth matter, neutrinos are, for coherent forward scattering in Earth matter, sensitive to 
the combination 

where e^^ = e^^ + e^^ . This is because the neutrino beams are only sensitive to the vector 
component. Further on, we will discuss how well one can test this combination. For the 
bounds on interactions for individual fermions, see Table 8 in Ref. [2]. 

2.1 Neutrino propagation Hamiltonian and considered NSI parameters 

Interactions of the type in Eq. ([T]) add an extra effective matter effect potential to the 
neutrino propagation Hamiltonian, which then reads 

'flee 
H = -^lU \ Am2 I f/^ + I I + 
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Here, ace is the usual matter effect term defined as ace = '^^^/^EGpNe (with N^ the electron 
number density in Earth matter), and the first line corresponds to the usual Hamiltonian 
in Earth matter. This equation already implies that the energy and baseline dependence of 
the non-standard effects will be similar to the standard matter effects, i.e., long baselines 
and high neutrino energies are important. For antineutrinos, the matter potential in Eq. ([3]) 
and all complex phases change sign, i.e., ace -^ "^cc, U -^ U*, and e™^ — > (e^^)*- Note 
that from the hermiticity of the Hamiltonian, e^^, are real numbers, while the e^a's can be 
complex for a ^ (3. As far as purely phenomenological bounds are concerned, |e^| and |ej]|j| 
are already very well constrained (see, e.g., Table 8 in Ref. [2]). In fact, we will show at 
the end of this study, that the bounds obtainable from the neutrino factory are comparable 
to the current bounds, which means that the neutrino factory is probably not the best 
experiment for their measurement. The interaction described by e™ is not per se interesting 
for us, since it will be intimately correlated with the matter density. We will discuss it in 
Sec. El 



^The coupling to the neutrino fields is left-handed because a right-handed coupling would be either 
helicity-suppressed, or only present in higher order corrections (there have to be at least two non-standard 
vertices in the amplitude to produce and absorb the right-handed neutrino). 



Because of the strong bounds on |e^| and |e"l |, and the straightforward relationship between 
e^ and the matter density precision measurement, we will focus on e™, e™ , and e^^ in the 
main line of this study. Note that the above mentioned bounds are purely phenomenological, 
and there are no convincing theoretical arguments yet why these non-standard effects should 
be large. Hence we focus on further constraints beyond the current limits in this study, but 
we do not discuss a possible discovery of non-standard effects, and only marginally touch 
possible effects on the determination of the standard oscillation parameters. Note that 
similar non-standard effects can be present in the neutrino production or detection. We do 
not consider these effects, which has the advantage that we do not have to simulate the near 
detector explicitelyo 

2.2 Measuring e^ in the golden and silver appearance channels 

Let us now first of all focus on e^, which can be best measured in the golden z/g -^ Ufj, and 
silver Ue -^ Vt appearance channels (see, e.g., Ref. [46]). The interference term induced by 
e^ in the silver channel P^t = Pv^^ur can be illustrated as 
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where "SO" stands for "standard oscillations". As we will see below, this structure is 
recovered in the full expression of the oscillation probability. In our discussion, we will use 
the following abbreviations for the spectral terms, i.e., the terms containing energy and/or 
baseline information: 

A - ^^. (5) 

J,,, ^ sin[(l - i)A] 

1-i ' ^^ 

a/2 
T^"" = sin(AA) = sin(±— Gi^iVeL) . (8) 

Here A corresponds to the vacuum oscillation phase, A to the effective matter potential 
with A ^ 1 at the matter resonance, J^^'^'^ to a term maximal at the matter resonance, 
and JF^^ to a term which is vanishing at the magic baseline L ~ 7 500 km [8,64]. In the 
definitions of A and JF^^, the upper signs are for neutrinos, and the lower ones for antineu- 
trinos. Contributions proportional to different products of these terms can, in principle, 
be disentangled by the use of a wide beam spectrum and different basehnes. The standard 



■^There is not yet any near detector specification in the IDS-NF baseline setup. As soon as such a 
specification is available, it may make sense to discuss production and detection effects as well. 



oscillation probability for Pe^ and Per is, to second order in a = /\m\^/ Am\^ ~ 0.03 and 
sin 26^13, given by (see, e.g., Ref. [65]) 
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sin 6ij and c,: 



^jj — cos 6ij. Note that the upper row/signs are for Pg^, and the lower 
row/signs for P^t. The different terms, can in principle, be disentangled by their spectral 
dependencies. For example, for jF'^^ —>■ (magic baseline) only the first term survives, 
which allows for a clean measurement of sin^ 26'i3 and the mass hierarchy. The relative 
amplitude of the different terms is given by the size of sin^ 26'i3 compared to a^ ~ 0.001: For 
sin^ 2^13 ^ a^, the first term dominates, for sin^ 26'i3 ~ a^, all terms including the middle 
(CP-terms) are large, and for sin^ 26*13 ^ «^) the last (solar) term dominates. Introducing 
non-standard effects by e^, we have to second order in a, sin 2^13, and |e 
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For antineutrinos, A changes sign, and all phases are inverted, i.e., the corresponding sin- 
terms change signs. That means that the 2nd and 5th terms in Eq. (TTO!) are the CP-odd 
terms describing intrinsic non-standard CP violation. Note that there can be CP violation 



even for sin^ 2^13 = 0, which is then induced by the 5th term. As we can read off from 
Eq. (jl]), there are terms proportional to |e^| sin 26*13, terms proportional to |e™|a, and 
terms proportional to |e^p, which dominate depending on the relative size of sin 2^13, a, 
and |e™|. For example, for sin^ 2^13 = 0, only the last six terms survive. If in addition 
|e™| :^ a ~ 0.03, the last three terms dominate, which are quadratic in |e™|. 



For the current best-fit value 623 = tc/4, P^jf and P^t differ only by the signs as given in 
the 2nd, 3rd, 5th, 6th, and 9th terms (and the 2nd and 3rd term in Eq. ([9])). If there is 
information from many different baselines and energies, the different dependencies on the 
spectral terms Eqs. ([S]) to ([H]) can be used to disentangle all terms in Eq. ( TTU]) except for the 
2nd and 5th (or 3rd and 6th) terms. Depending on the relative size of sin 2^13 and a, either 
of these two terms may dominate, or both terms might be of similar magnitude. Note, 
however, that in certain limits, Eq. (ITOl) is very different for the golden and silver channels. 
For example, let us consider the situation at peak energies of the spectrum, and at the first 
oscillation maximum, which occurs typically at a baseline around L ~ 3000 to 4000 km. In 
this case, A A ^ tt/2 and A ^ 1. It is easy to see from Eq. i^ to Eq. ([H]), and from Eq. 1^ 
that this leads to all standard oscillation terms being oc A^ oc l/E"^, so that they cannot 
be disentangled from each other. Of the NSI terms in Eq. (TTOj) . only those which are either 
constant in energy, or proportional to 1/E, may be separated from the standard terms. 
From the signs in Eq. (fTOl) . we find that, for the golden channel, all relevant NSI terms 
cancel, while for the silver channel, they interfere constructively. Consequently, the golden 
channel detector at the short neutrino factory baseline will not be able to provide a good 
sensitivity to e™ , while a better performance is expected for a silver channel detector at the 
same baseline. Our argument also shows that a high neutrino energy is advantageous for a 
measurement of NSI at AA ~ n/2 and A <C 1, because it reduces the standard oscillation 
background. 

At the magic baseline, JF'^^ — > [8], we obtain 



pNSI \ / e2 \ 



+ A\eZ\ sin20i3C23 ( 'f ) cos(5cp + C) ^ (•^'''Imb)' 

\ '^23 / 

+ 4|eg^| C23 I 2 1^1-^ ImbJ • (iij 

Here J^^'^'^Imb is J^^"^^ in the magic baseline limit A A -^ n. This formula is exactly the 
same as in Ref. [51] for Pg^ if -^^'^'^Imb is trigonometrically expandedo It has a number of 
interesting implications. First of all, there are much less correlations than in Eq. flTOj) . which 
means that the magic baseline will crucially contribute to the NSI sensitivity. However, 
compared to the SO case, the NSI case is not completely correlation-free at the magic 
baseline - even the phase 0™ appears in the formula. Second, for maximal atmospheric 
mixing, we have Pe/x = Per, which means that there is no difference between the golden and 



*We keep, however, J-"^°''|mb in the formula, because we can even correctly reproduce the resonance limit 
i ^ 1. In Ref. [51], Pe^ -> 00 for i ^ 1. 



silver channels. Therefore, there will be no physics case for the silver channel at the magic 
baseline because of the much lower event rate. And third, since all terms are proportional 
to {J-'^^^)'^, the second term, which is proportional to A oc E, and the third term, which is 
proportional to A^ oc E"^, become relatively enhanced for high energies. This means that 
high neutrino energies are very important to constrain NSI. From Eq. flTTl) . we can already 
estimate that the |e^p sensitivity should quantitatively be comparable to the sin^ 26'i3 
sensitivity, i.e., if the sin^ 26'i3 sensitivity is about 10~^, we obtain a |e^| sensitivity of 
about 0.003 if correlations and degeneracies can be sufficiently resolved. 



2.3 Measuring e^ and e"^ in the disappearance channel 

As we shall quantitatively discuss later, the disappearance channel P^^ at the neutrino 
factory is the dominant source for the ej^ and e^^ sensitivities (see, e.g., Ref. [47]). Here 
we follow Ref. [44] to describe these effects in the two flavor limit. The approximation 
corresponds to the u^-Ur system with 6*13 -^ 0. For e"^^, we have 



where the PMNS matrix [/ is a 2 x 2 mixing matrix with the mixing angle corresponding 
to 6*23. From this expression, we can read off the fact that e!^ plays the same roll as ej^ 
does [47]. In this case, we can describe the shift in the mass squared difference and mixing 
angle by a parameter mapping: 



Arhl-^ = Am^^ Wsin^ 2^23 + (ierr + cos 2^23) , (13) 

sm 26'23 = 7^ (14) 

sin2 2e23+ fie™ + cos 2^23 

In the maximal mixing limit ^23 -^ 7r/4, they are reduced to 



Ami, ^ Am2,^l+(ie-)', (15) 

sin2 2^23 = r^. (16) 



The lowest order of this shift comes from (9{(e^)^}, which means that it does not appear in 
the analytic expressions in Refs. [52,53]. The NSI effect is proportional to A'^ oc E"^. In low 
energy experiments such as T2K, this effect is not important. On the other hand, in high 
energy experiments, such as neutrino factories, this will affect the oscillation probability 
signiflcantly. In addition, note that there can, in principle, be resonant effects for strong 
deviations from maximal mixings. For |e™ | = 0{1) (which might be, however, unrealistically 
large [29]) and ^23 on the edge of the current 3a allowed range, one flnds from Eq. (TT^ that 
one can have resonance energies as high as about 2.5 GeV, which is slightly above the 
currently considered detection threshold. 



For e"^, which is also strongly present in the disappearance channel, the parameter mapping 
is slightly more complicated because e™ can have a complex phase 0!^: 

Am^^ = Am2,^(^2i|e-|cos0™ +sin2^23)'+ (2i|e™ | sin</);f^)' + cos^ 2^23, (17) 



(2i|e;^,| cos0;7, + sin 2^23)' + (2i|e;^,| sin^;^,' 



fl^ 
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(19) 



2i|e;f^|cos0™+ sin 2^23) + (2116™ | sin 0™) +cos22^; 
For maximal mixing ^23 -^ 7i"/4, we obtain sin^ 2^23 -^ 1 and 

Ami, -^ Am^.^l + 4i|e™ | cos0™ + (2i|e™ |)' . 

We can see that the mass squared difference receives modifications already at first order in 
e™. , while the mixing angle remains maximal to all orders. Since we will marginalize over 
the phase of the NSI parameter, the visible effect comes from the second order term. 

Numerically, the sensitivity to |e]^| and |e™ | will be limited by the precision of Amli 
(provided that all other correlations can be resolved). In Ref. [9], the la precision of Amli 
has been found to be 0.2% at the very long baseline. Eq. (TTSll thus implies 0.5 (A|e!^|)^ ~ 
0.2% at the la sensitivity limit for |e^|. At the upper end of the neutrino spectrum 
(E^ ^ E^ = 25GeV, A ~ 3), this leads to le"^^] ~ 0.02. From Eq. ([19]), we obtain a much 
better sensitivity for real e™ , i.e., cos0"^ = ±1: In this case, the sensitivity is linear in 
|ej^|, and given by 2yl|e™ | ~ 0.2% at the sensitivity limit, or \e^^\ ~ 3 ■ 10""^. If, however, 
0™ can take any value, it can also assume 0™ = ±7r/2, and we are back in the quadratic 
regime such as for |e^|. In fact, one can even have cancellation of the two terms in Eq. flT9|) . 
which means that we expect a sensitivity worse than for |e™ |. 

3 Performance indicators and simulation details 

In the previous section, we have motivated why we only consider small non-standard effects. 
As performance indicator, we use the "|e^fl| sensitivity", which corresponds to the exclusion 
limit which is obtained if the true value (simulated value) vanishes. In principle, we follow 
the same definition as for the sin^ 2^13 sensitivity, i.e., we define the |e^^| sensitivity as the 
largest fit |e^^| which fits the true |ej^^| = 0. Note that e™^ can be complex for a ^ (3, 
which means that the (fit) phase 0^^ has to be marginalized over, whereas the true phase 
is irrelevant because of the true |e™o| = 0. In addition, all standard oscillation parameters 
are marginalized over. In our simulations, we find the main correlation leading to technical 
difficulties is the correlation among 0J^^, (5cp, and sin^2^i3. Therefore, we pre-scan this set 
of parameters in many cases to find the position of the global minimum. In some cases, 
we will consider also correlations among different ej^^'s in order to compare our results to 
earlier works. For the same reason, we will sometimes also neglect the phases 0^^ even for 
a ^ (3. For the sake of simplicity, we do not include the sgn(A?Ti|^) degeneracy for the non- 
standard sensitivities [11]. In addition, we do not consider degeneracies with unrealistically 



large |e™^| > 1 in some cases, because these degeneracies would appear above the current 
bounds. Note that our |e^^| sensitivity is expected to be similar to a conservative case 
discovery limit, i.e., depending on the phases, the discovery may be possible for smaller 
|ej^^| than the sensitivity hmit (c/., Refs. [50,52]). 

The experimental scenario we consider is the IDS-NF 1.0 setup from Ref. [63], which is the 
current standard setup for the "International design study of the neutrino factory" (IDS- 
NF). Within the "International scoping study of a future neutrino factory and super-beam 
facility" [2,9], this setup has been optimized for the measurement of sin^ 2^13, the neutrino 
mass hierarchy, and leptonic CP violation in the case of standard oscillations. In short, this 
setup uses two baselines at about 4 000 km and 7500 km with two (identical) magnetized 
iron neutrino detectors (MIND) with a fiducial mass of 50 kt each. In addition, a 10 kt 
emulsion cloud chamber (ECC) for z/^ detection is placed at the short baseline. For each 
baseline, a total of 2.5 ■ 10^^ useful muon decays plus 2.5 ■ 10^^ useful antimuon decays 
in the straight of the corresponding storage ring is used, which could be achieved by ten 
years of operation with 2.5 ■ 10^'^ useful muon decays per baseline, year, and polarity. The 
muon energy i?^ is assumed to be 25 GeV, which is sufficient for a detector with a low 
enough detection threshold [9]. The detector and systematics specifications can be found 
in Refs. [7,63]. Note that there is not yet any near detector specification. We do not 
simulate the near detector explicitely, because we do not discuss non-standard production 
or detection effects such as in Ref. [52]. In addition, we do not require charge identification 
in the disappearance channel, which means that we have to add the u^ and z/^ event rates. It 
has been demonstrated in Ref. [9] that the better efficiencies (and better energy threshold) 
lead to a better performance in that case. In summary, the following oscillation channels 
are included: 

u^ at 4 000 km (z/^ appearance) 
z/^ at 4 000 km (z/^ appearance) 
3. z/g -^ z/^ at 7500 km (z/^ appearance) 
z/^ at 7500 km (z/^ appearance) 
J^e ^' i^fj. + ^^JL at 4 000 km (z/^ disappearance) 
t'e ^ i^/i + J-'fi at 4 000 km (z/^ disappearance) 
t'e ^ i^/x + J^fi at 7500 km (z/^ disappearance) 

8. i>/x + z/g — !► i>yx + z/^ at 7500 km (z/^ disappearance) 

9. z/g -^ Vt- at 4 000 km (z/^ appearance) 

In the following, we will refer to golden channels 1 to 4 as Golden [5], to channels 5 to 8 
as the disappearance channels, and to the silver channel 9 as Silver [6]. Note that in the 
limit of small sin^2^i3, channels 5 to 8 can be approximated by u^ -^ u^ and u^ -^ u^, 
respectively, which we have discussed in the phenomenology section. 
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Compared to Ref. [63], we study several modifications in order to discuss the neutrino 
factory optimization for non-standard interactions. In some cases, we will vary the muon 
energy or baseline(s), sometimes even for Silver separately. In addition, we will discuss a 
potentially improved silver channel detector Silver*, which uses five times the signal and 
three times the background of Silver in order to implement the hadronic decay channels of 
the r as well [9]. In a part of the study, we will not include Silver at all. 

All simulations are performed using the GLoBES software [66,67]. The experiment de- 
scription is based on Refs. [9, 12] updated with the numbers from Ref. [63]. For the 
true oscillation parameters, we use sin^ 6*12 = 0.3, sin^ ^23 = 0.5, Amli = 7.9 ■ lO^^eV^, 
Amg^ = 2.6 • 10^'^ eV^, and a normal mass hierarchy unless stated otherwise. For the true 
sin^ 2^13 and 6cp, we choose certain benchmark points, but we will see that there is relatively 
little dependence on their true values in most cases. For 612 and Am2i, we assume external 
measurement precisions of 10% each, whereas we do not impose any external constraints on 
the leading atmospheric parameters. The used true values and their errors are motivated 
by the current best-fit values and their errors, see, e.g., Refs. [1,68]. For the matter density, 
we use the PREM profile (Preliminary Reference Earth Model) with a normalization uncer- 
tainty of 5% [69,70]. For neutrino trajectories which do not cross the core of the Earth, we 
approximate the PREM profile by a single layer of constant density, while for core-crossing 
neutrinos, we use a mantle-core-mantle profile with three layers. The densities within the 
respective layers are computed by averaging the full PREM profile along the neutrino tra- 
jectory. Note that the 5% matter density uncertainty is assumed to be correlated among 
different channels operated at the same baseline, and uncorrelated between different base- 
lines (unless we vary the baseline of one channel independently; in that case, it is always 
uncorrelated) . 

4 Impact of different channels 

In this section, we discuss the impact of different oscillation channels, and we study the 
optimization of the silver channel. We know from Ref. [51] that the combination of two 
baselines, one with about 3 000 km and the other with about 7000 km, turns out to be 
very useful to resolve correlations between the standard and non-standard parameters, and 
among different non-standard parameters. However, disappearance information was not 
taken into account in the analysis of Ref. [51], and the off-diagonal e's were assumed to 
be real. On the other hand, it has been demonstrated in Ref. [46] that the silver channel 
probability at 3 000 km significantly depends on the non-standard effects, especially e™. 
Therefore, we focus on three major questions in this section: 

1. Which oscillation channels dominate the measurements for which non-standard quan- 
tities? 

2. If one has already a two-baseline setup, such as the IDS-NF setup, does one still need 
the silver channel? 

3. Is the silver channel location at the shorter of the two golden baselines really the 
optimal choice? 
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Figure 1: Contribution of different channels to the e^-e™ sensitivity. The first row corresponds to 
the golden channel only. The different columns correspond to the two different baselines 4 000 km and 
7 500 km, as well as their combination. In the second row, we in addition add the disappearance channel 
(left), introduce complex e™ (middle), and finally add the Silver* channel (right). In the upper row, we only 
marginalize over sin 1d\j, and (5cP) whereas in the lower row, we marginalize over all oscillation parameters. 
Note that e™ is assumed to be real in the first four panels, and complex in the last two. In this figure, a 
true ^cp = 37r/2 and sin^ 2^13 = 0.001 have been assumed. In addition, E'^ — 50GeV has been chosen for 
comparison to Ref. [51]. The contours correspond to the Icr, 2cr, and 3o" confidence level for 2 d.o.f. 

These questions can only be quantitatively and reliably answered using a full simulation. 

To compare our results to Ref. [51], let us first of all assume all e^^ to be real. In addition, 
we study simultaneous constraints for two non-standard parameters to illustrate the impact 
of different channels. In order to compare to Ref. [51], we choose an example in the 6™-e™ 
plane, where correlations are particularly severe, i.e., sin^ 2^5^3"'' = 0.001 and S^^ = 37r/2. 
In Fig. [T] we show the allowed sensitivity region in the upper row for L = 4 000 km (left 
panel), L = 7500 km (middle panel), and the combination of the two baselines (right panel). 
In these panels, we have only marginalized over sin^ 26*13 and ^cp, the true e™ = e™ have 
been assumed to vanish, and we have chosen E^ = 50 GeV for the whole figureO Our results 
reproduce Ref. [51] very well, even though the baselines are slightly changed to match the 
IDS-NF baseline setup. Note that the magic baseline fit is not completely correlation-free. 



^We will discuss the impact of the muon energy in the next section. 
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as it is obvious from Eq. (TTTI) . which means that there is no clean measurement of e™ 
at the magic basehne. In addition, there are still some problems with correlations in the 
combination of the two baselines. Therefore, one may suspect that the silver channel could 
help to resolve these. 

We have tested that this problem becomes even worse if one marginalizes over the leading 
solar and atmospheric oscillation parameters as well. In this case, the silver channel indeed 
helps to resolve the correlations. However, this conclusion does not hold anymore if one in 
addition adds the disappearance channel, as we have done in the lower left panel of Fig. [1] It 
helps to measure the atmospheric oscillation parameters, and it severely constrains e^ (or, 
as we have tested, e"^). Furthermore, there is almost no correlation remaining between e"^^ 
and e^. From the comparison with the upper right panel we learn that the golden channel 
indeed has the best e^ sensitivity, whereas the disappearance channel has the best e^^ and 
e^ sensitivity. Therefore, the analytical formulas presented in Sec. |2] are really the ones 
applicable to the most sensitive channels at the neutrino factory. 

Now one can argue that once even more parameters are added, the information from the 
silver channel needs to contribute at some point. In addition, we have not included the 
complex phase 0^ yet. Therefore, we show in Fig. [H lower middle panel, the full complex 
case with the additional parameter 0^ marginalized over. Note that the scale on the axes 
has changed, and that e^^ is real by definition. There is no large quantitative change 
compared to the lower left panel. However, if the Silver* channel is in addition used at the 
4 000 km baseline, as we illustrate in the lower right panel, |e^| can be somewhat better 
constrained, whereas there is almost no effect for e™ (or e"^^). Though this improvement is 
not insignificant (a 30% effect), we will study in Sec. [5] how it quantitatively depends on 
the muon energy and silver channel implementation. We have checked that it cannot be 
achieved by a mere up-scaling of the Golden detectors, as one may naively expect, even if 
the detector masses are increased by 50 kt each. Therefore, we have identified a synergy 
in the sense of Ref. [71] here. We have checked that even if one includes in addition e"^ 
to be marginalized over, there are no significant qualitative changes to this picture, i.e., no 
additional correlations to be resolved. Therefore, we conclude that the Silver* channel at the 
4 000 km baseline is not a key component to push the non-standard parameter measurements 
by an order of magnitude, but it may help to improve the |e^| sensitivity somewhat. One 
reason are the relatively low event rates even for 6cp = 37r/2, where the silver rate becomes 
largest in the CP-odd term (second term in Eq. (jH])): in total about 47 events for z/,- 
appearance (Silver*) or 9 events for Ur appearance (Silver), compared to about 323 events 
for i>n appearance and 6 million events for u^ disappearance at the 4 000 km baseline (i?^ = 
50 GeV, sin^ 26'i3 = 0.001, 6cp = 37r/2, normal hierarchy, no NSI). On the other hand, we 
know from the upper left and middle panels in Fig. [T] that the magic baseline significantly 
contributes to the e^ measurement for the golden channel in an orthogonal way, and there 
is still a substantial number of events at this baseline (126 events for the above benchmark 
point). Therefore, we have demonstrated that the silver channel at the short baseline is not 
mandatory for the e^ sensitivity if the golden channel at the magic baseline is used. In 
addition, at the magic baseline, the golden and silver appearance channels are equivalent, 
as we have found in Eq. flTTl) . Hence, we do not expect the silver channel to be useful at 
the magic baseline either. 



12 



Even if the silver channel does not help at 4 000 km or 7 500 km, what about (hypothet- 
ically) placing the ECC at a third baseline in combination with two Golden detectors at 
4 000 km and 7500 km? We discuss this question for |e^| in Fig. [2] for two (representative) 
sets of sin^ 29^^^'^ and (^c™*^. The different curves correspond to the different muon energies 
25 GeV, 50 GeV, and 100 GeV. Note that in this case, the matter density is assumed to 
be uncorrelated among all three baselines. Indeed one can read off these figures that the 
optimal performance is obtained at about 3 000 km to 4 000 km, where the vertical lines 
correspond to our standard choice. This means that our setup is perfectly optimized for the 
silver channel. In addition, one can read off this figure that for E^ ^ 50 GeV hardly any 
effect is visible. We will discuss the muon energy dependence in detail in the next section. 

Note that we will not consider two NSI parameter correlations for the rest of this study 
anymore, since we have found that there is hardly any correlation remaining in the e™-e^ 
plane (c/.. Fig. [H lower middle and right panels). This means that the sensitivities can as 
well be studied separately^ Of course, there is a straightforward correlation in the ej^-e^ 
plane, which can be analytically understood from the full two-parameter mapping [44] . For 
example, for maximal atmospheric mixing and real e"^, we have (c/.. Sec. 12.31 and Ref. [44]) 



/A 



~ 2 

m 



(2^6- + 1)^ + (Ae- )^ = (^^J = const. , (20) 

which just corresponds to the normal form of an ellipse centered at e™ = — 1/(2A) and 
e^^ = 0. We therefore do not consider this correlation anymore. The e™-e™ plane, on 
the other hand, is similar to the e™ -e^ plane, as it is obvious from Sec. I2.3[ In principle, 
correlations between e™ and the other parameters could be interesting, but this would be 
out of the main line of this study. Note that the correlation with e^ is intimately connected 
to the matter density uncertainty, which we have taken into account. This means that our 
simulation using a 5% matter density uncertainty is, apart from the fact that the matter 
density is assumed to be uncorrelated between the two baselines, equivalent to a simulation 
with a precisely known matter density profile and an external bound |e^| < 0.05 (Icr). We 
will further comment on the relationship to the matter density uncertainty in Sec. [51 

5 Optimal muon energy 

From Fig. [2], we have learned that the muon energy has some impact on the NSI performance. 
This can, for example, be seen at the relative dominance of the NSI terms for higher neutrino 
energies in Eq. fITT]) . Therefore, we discuss in this section the dependence of the NSI 
parameter sensitivities on the muon energy, with and without silver channel. Note that 
for the standard oscillation parameters, the muon energy has only a minor impact if the 
detection threshold in the Golden detectors is low enough and E^ > 20 GeV [9]. Therefore, 
a lower muon energy Efj_ = 25 GeV was chosen for the IDS-NF baseline setup than the 



^In fact, for very large sin 29i3 close to the current bound, there is some correlation in the e^'^-e™ plane 
remaining, which partly comes from the matter density uncertainty. Since for sin 20i3 > 0.01, a neutrino 
factory would probably look different from the current IDS-NF baseline setup (i.e., have a short baseline 
and a lower muon energy) [72-74] , we do not discuss this case anymore. 
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Figure 2: Sensitivity to |e^| (3cr) as a function of the Silver* baseline for three detector setups: two 
MIND detectors at 4 000 km and 7 500 km, and one Silver* detector at the specified basehne on the horizontal 
axes. The different curves correspond to different muon energies as given in the plots. The different panels 
show the result for different (representative) true values of sin 2^13 and (5cp, as given in the captions. The 
IDS-NF iron detector short baseline is marked by the vertical lines. 

originally anticipated E^ = 50 GeV. For a possibly different detector technology, also a 
low energy neutrino factory with E^ r^ b GeV has been discussed in the literature [72-74]. 
Such an experiment may be useful for large sin^ 26'i3, but the typical setups involve only one 
baseline, and will therefore not be considered in this work. Again, there are three relevant 
questions for this section: 

1. Is the muon energy of 25 GeV sufficient for the NSI sensitivities, or should one go to 
a higher E^ for the IDS-NF baseline? 

2. What are the prospects to improve the current NSI bounds for a considerably lower 
muon energy? 

3. How important is the silver channel as a function of the muon energy? 

In order to address these questions, we present in Fig. [3] the sensitivity to |e™|, |e™ |, and 
|e™ I (3o") as a function of the muon energy for two different sets of true values. The 
different curves correspond to different |e^^| and different detector setups: two magnetized 
iron detectors at 4 000 km and 7500 km, and one optional Silver or Silver* detector at the 

and 5cp°' which, however, 
\, however, we hardly find 



shorter baseline. For |e™ |, we find some dependence on sin^ 29^^^^ 



HT 



does not lead to qualitatively different conclusions. For e 
any dependence on sin^ 26^^^^ and (Jqp'^. 

Neglecting the silver channel for the moment, we find in all cases in Fig. |3]a strong depletion 
of the sensitivity for E^ < 20 GeV. For higher muon energies, however, there is no significant 



14 




40 60 

E, [GeV] 



GLoBES 2008 
80 100 




40 60 

E^ [GeV] 



GLoBES 2008 
80 100 



Figure 3: Sensitivity to |e™|, |e!J^^|, and |e"i^| as a function of the muon energy for different detector 
setups (3cr): two magnetized iron detectors at 4 000km and 7 500km, and one optional Silver or Silver* 
detector at the sliorter baseline. The two panels show the result for different (representative) true values of 
sin^ 2^13 and (5cp, as given in the plots. Here a normal mass hierarchy is assumed. The IDS-NF standard 
muon energy (25 GeV) is marked by vertical lines. 



gain anymore. Ttie reason is that the energy range with the strongest matter effects is 
sufficiently covered, while for higher muon energies, the event rates at the lower end of 
the spectrum decrease somewhat. In fact, the |e^| sensitivity even becomes worse for 
E^ > 50 GeV. This means that E'^ = 25 GeV is indeed sufficient for the NSI sensitivities. 
For E"^ ~ 5 GeV, however, the current bounds could only be improved for |e™ | and |e^| by 
a factor of a few, and not at all for |e"^|. Note that also the standard oscillation parameter 
measurements are significantly affected for E^ <^ 20 GeV [9]. 

For the silver channel, we do not find any impact for the |ej^| and |e^| sensitivities what- 
ever the chosen muon energy is. For the |e^| sensitivity, however, it slightly improves the 
sensitivity for Efj, > 25 GeV, depending on the true parameter values. For E^ = 25 GeV, the 
silver channel hardly contributes. This means that, at least for the NSI, the current choice 
of E^ = 25 GeV for the IDS-NF baseline setup is in contradiction with the Silver detector at 
4 000 km. Either the muon energy needs to be increased to make the silver channel valuable, 
or there is no physics case for the silver detector. Since the impact of the silver channel is, 
even for E^^ = 50 GeV, not very large, we prefer to choose E'^ = 25 GeV in the following 
in order to be consistent with the IDS-NF baseline setup. However, we will not include the 
Silver detector, because we have not found any significant physics contribution. 

As far as the |e™| sensitivity is concerned, let us first of all make some numerical estimates. 
From Eq. ([3]), we know that ace + '^cc ^^ enters the Hamiltonian in the presence of e™. 
If, in addition, a matter density shift x is considered, i.e., ace -^ ctcc + x, the matter 
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Figure 4: Sensitivity to |e™| as a function of the muon energy for two magnetized iron detectors at 
4 000 km and 7 500 km {3a). The left panel is computed for a matter density uncertainty dace — 5%, and 
the right panel for a matter density uncertainty Sacc = 1% (!'''). Here different values for sin^ 26\'^° and 
(5qp° have been used, as shown in the plot legend. The IDS-NF standard muon energy (25 GeV) is marked 
by vertical lines. 



potential ace is (to leading order) shifted hj P = x + acc^Te- ^'^^ ^Te — 0' '^•^•' without 
NSI, P = X is exactly the deviation of the matter density from the reference matter density 
profile. Note that we impose a prior on the shift x, which means that |a;| < Sacc is limited, 
where 6acc is the matter density uncertainty (Icr). This uncertainty comes, for example, 
from the limited precision of seismic wave experiment. The measurement of P, i.e., the 
matter density precision measurement, has been studied in the literature in Refs. [75-77]. 
Therefore, one can, in principle, estimate the e™ sensitivity from the precision of P and 
Sacc as the Gaussian average for a single baseline experiment. For two baselines, e™ will 
be correlated between the baselines, whereas x will be not. Therefore, a slightly better e™ 
sensitivity might be expected in practice. 

We show in Fig. H] the sensitivity to |e™| as a function of the muon energy for two magnetized 
iron detectors at 4 000 km and 7 500 km {3a). The left panel is computed for a matter density 
uncertainty Sacc = 5%, and the right panel for a matter density uncertainty Sacc = 1% 
(la). Here different values for sin^2^5^3"® and S^p"^ have been used, as shown in the plot 
legend. For very high E^, the effective precision on P will be higher than the matter density 
uncertainty Sacc, which means that the |e™| sensitivity will be asymptotically limited by 
Sacc- Therefore, in the left panel of Fig. HJ the sensitivities are roughly limited by 3 ■ 
0.05/v2 ^ 0.11, in the right panel by 3 ■ O.Ol/y^ ~ 0.02, where the factor y^ comes from 
the fact that two independent matter density priors are added. For very low Efj_, the precision 
of P will be much weaker than Sacc in all cases, which means that the sensitivity is limited 
by the precision of P. Therefore, the curves in both panels are very similar for small muon 
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energies irrespective of the matter density uncertainty. The dependence on sin^ 29f^^ is 
similar to that of the measurement of P without NSI, see Ref. [76]. Note that there is hardly 
any dependence on 6qp^ for very small or very large sin^ 29f^^. For very large sin^ 2^13, the 
precision of P is extremely good already for comparatively small E"^, which means that the 
asymptotic limit is quickly reached (c/., curves for sin^ 26*13 = 0.1). For very small sin^ 29]^^^, 
the solar term (fourth term in Eq. iQ) dominates the measurement, which means that the 
performance is poorer than in the large sin^ 26*13 limit |j For intermediate sin^ 26*^3"'^, the 
performance strongly depends on ^^p^, because Sep leads to non-trivial correlations (c/., 
curves for sin^2^j3"'^ = 0.001). Compared to the other sensitivities discussed in this section, 
there can be a strong gradient between a muon energy of 25 and 50 GeV especially if 6acc 
is sufficiently small, which, however, somewhat depends on the true parameters. Therefore, 
if one emphasizes the |e^| sensitivity, a higher muon energy might be important. Since 
the improvement would only be a factor of a few beyond the current bounds, this would 
probably not be the main argument for a higher muon energy. 

6 Baseline optimization: Standard versus non-standard physics 

In this section, we discuss the optimization of the baselines for a neutrino factory with 
two detectors, considering both standard oscillation physics and non-standard scenarios. 
The experimental setup is based on the IDS-NF 1.0 configuration, but omitting the silver 
channel. We will treat the two baselines Li and L2 as free parameters in this section, 
and compare their resulting optimal values to the ones suggested by the IDS-NF 1.0 setup, 
namely Li ~ 3 000 km to 5 000 km, and L2 ~ 7 000 km to 8 000 km. For definiteness, we 
define our benchmark setup by Li = 4 000 km and L2 = 7 500 km. This choice is based on 
previous works, in particular on the magic baseline argument in Ref. [8] (correlations and 
degeneracies disappear at a baseline around 7500 km), the single baseline optimization in 
Ref. [9] (leading to the conclusion that one needs two baselines to optimally measure all 
standard performance indicators), and the optimization of the longer baseline L2 with the 
shorter one fixed at Li = 4 000 km in Ref. [77]. A simultaneous variation of both baselines 
has so far only been considered in Ref. [8] for the sin^ 2^13 sensitivity reach. 

The main questions relevant for this section therefore are: 

1. Is the IDS-NF baseline setup still optimal for all standard oscillation performance 
indicators if both baselines are varied simultaneously? 

2. Is the standard optimization robust if there are NSI? 

3. What would the optimal baselines be for non-standard interactions? 

4. Is the non-standard optimization consistent with the standard one? 



^In fact, the fourth term in Eq. ([9]) is CP-invariant, which means that it is also invariant under e™ = 
-^ e™ = — 2 (which corresponds to a sign flip of the matter density profile). We have not included this 
additional degeneracy. 
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Optimization for standard oscillation performance indicators 

In order to discuss the optimization for the standard performance indicators, we show in 
Fig. O the two-baseline optimization of a neutrino factory with two Golden detectors. In 
this figure, optimal performance means optimal reach in sin^ 26'i3. The upper left panel 
shows the region with optimal sensitivity to sin^ 2^13, the upper right panel shows the 
region with optimal sensitivity to the normal mass hierarchy, the lower left panel shows 
the region with optimal sensitivity to maximal CP violation, and the lower right panel 
shows the intersection of the three regions. The contours have been chosen such that all 
regions are of similar size (see values in plots). All performance indicators are defined at 
the 5a confidence level in order to include all degeneracies (even if they occur only at a 
relatively high x^)- For the sensitivity to the mass hierarchy and to CP violation, we use 
Sqp^ = 37r/2, since for this value, degeneracies have a strong impact, so that it corresponds 
to a conservative assumption. The circles mark our IDS-NF benchmark setup, whereas the 
colored (gray) diamonds mark the optimum baselines for each case. The dotted curves and 
the black diamonds indicate the optimization in a non-standard scenario, which will be 
discussed below. 

Let us, however, first focus on the standard optimization only. For the sin^ 2^13 sensitivity 
(upper left panel), we recover the shape from Ref. [8], Fig. 1. In this case, almost any 
combination of baselines, with one of them being magic, leads to a good performance. The 
same conclusion can be obtained for the mass hierarchy sensitivity (upper right panel), where 
the the optimal baseline somewhat varies with ^^p'' (c/., Ref. [9]). Note that for sin^ 26*13 as 
well as for the mass hierarchy, the long baseline is a prerequisite. For CP violation (lower left 
panel), however, one of the baselines has to be short, i.e., between 2 000km and 6 000km, 
while a long baseline is required to resolve the degeneracies (for ^^p*^ = 7r/2 one might not 
need that [9]). Interestingly, no points on the diagonal are within the optimal region, which 
means that it is not sufficient to use only one baseline for this set of parameters. Finally, 
the lower right panel shows the intersection of the other regions. It clearly demonstrates 
that our standard choice, denoted by the circle, is well within the optimal region for all 
performance indicators. 

As the next step, assume that the non-standard effects are taken into account in the fft, 
which will obviously spoil the standard oscillation parameter sensitivities of a neutrino 
factory (see, e.g., Refs. [34,41]). This means that we have to marginalize over these effects 
as well. Let us first of all focus on e^: The dotted curves in Fig. O have been obtained from 
a fit where also e^ (its absolute value and phase) has been marginalized over. As one can 
read off the contours, the absolute performances for all the standard oscillation parameters 
become worse. However, the optimization does not change, as it can be read off from the 
lower right panel. Therefore, our two-baseline optimized setup is very robust even with 
respect to non-standard e^. We have also checked the other e™^'s for this optimization. 
While e^, e^^, and e™ hardly have any effect on the appearance channel at all (see also 
analytical formulas in Ref. [52]), e™ has a similar qualitative effect as e^. However, if 
the stringent existing bounds on |e^| are taken into account, the possible effects of this 
parameter become completely negligible. 
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Figure 5; Two-baseline optimization of a neutrino factory (with two Golden detectors) for the standard 
oscillation performance indicators. The upper left panel shows the (shaded) region where the sensitivity to 
sin^ 2^13 better than 10^^'^ i^'^): the upper right panel shows the (shaded) region where the sensitivity to the 
normal mass hierarchy (MH) is given for all sin^ 26*13 > 10~'^'® (5cr, S^p^ = 37r/2), the lower left panel shows 
the (shaded) region where the sensitivity to maximal CP violation (CPV) is given for all sin^ 26*13 > lO^"^'^ 
{5a, Sqp'^ = 3tt/2), and the lower right panel shows the intersection of the three regions as the shaded 
region. The dotted curves have been obtained from a fit including e™ marginalized (for the sin^ 20i3 ranges 
given in the plots). The diamonds show the setups with optimal sensitivities (colored/gray for the shaded 
contours, black for the dotted contours), whereas the circles correspond to the IDS-NF standard choices 
Li = 4 000 km and L2 = 7 500 km. 
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Optimization for the non-standard sensitivities 

We discuss the two-baseline optimization of a neutrino factory (with two Golden detectors) 
for non-standard interactions in Fig. O In this figure, the different rows represent the non- 
standard interaction parameters |e^|, |e™ |, and {e^^-l, respectively, whereas the different 
columns represent E^ = 25 GeV (left) and E^ = 50 GeV (right). Let us focus on E^ = 
25 GeV (left column) first, which is the IDS-NF choice. In these figures, the dashed curves 
show the optimum configurations for the standard optimization (shaded region from the 
lower right panel of Fig. \5^ for comparison. For the |e^| sensitivity (upper left panel), 
the IDS-NF baseline combination 4 000km plus 7 500km is close to optimal for both the 
standard (circles) and non-standard (diamonds) sensitivities. In this case, slightly longer 
rather than shorter baselines are preferred. Note that here the main contribution comes 
from the appearance channels. For the \e'^T-\ and |e^| sensitivities (middle and lower left 
panels), respectively, much longer baselines are preferred, such as core crossing baselines 
L > 10 700 km. The standard and non-standard optimizations do not coincide, which means 
that one would (hypothetically) need a third baseline. Here the main contribution comes 
from the disappearance channels, which tend to perform better at long baselines because 
more oscillation nodes can be resolved. 

In the right column of Fig. [6|, we increase E^j, to 50 GeV. While the absolute sensitivities 
at the standard choice 4 000 km plus 7500 km remain almost unaffected, better absolute 
sensitivities can be obtained for longer baselines. For example, for the |e™| sensitivity, 
the bound could be improved by about a factor of two if one went to a different (longer 
baseline) combination. The reason are the higher event rates for E^ = 50 GeV, which allow 
for better statistics at even longer baselines. Note, however, that the optimal region for e™ 
somewhat depends on the chosen sin^ 26*^3'^'^ and S^p"^, whereas the ones for ej^ and e^ are 
almost independent of these parameters. Since the qualitative discussion does not change, 
we decided to present the results only for one set of parameters. In addition, we have checked 
the optimization for the \e^\ and the |e™ | sensitivities. For the |e"^| sensitivity, there are 
hardly any qualitative and quantitative changes compared to the |e^| sensitivity (see also 
analytical discussion in Sec. 12.31) . For the e^ sensitivity, however, one baseline should be 
rather short 2 000 km ^ -^ ^ 4 000 km whereas the other can be long. This means that 
the standard IDS-NF baseline choice is close-to-optimal for this sensitivity. However, the 
optimal absolute sensitivity is about 0.005 {3a), i.e., not better than the current bounds, 
which means that this aspect is of little relevance. 

In summary, at least one very long baseline is an important prerequisite to put stronger 
bounds on the non-standard interactions. For Efj_ = 25 GeV, the standard and |e™| opti- 
mizations are consistent, while |e!^| and \e^^\ prefer one baseline to be as long as possible. 
However, even better absolute sensitivities could be achieved for longer baselines in combi- 
nation with a higher muon energy. 
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Figure 6: Two-baseline optimization of a neutrino factory (with two Golden detectors) for non-standard 
interactions. The different rows represent the sensitivities to the non-standard interaction parameters 
|e™|, \e"\.\, and |e"V|, respectively {3a), whereas the different columns represent Ef^ — 25GeV (left) and 
E^ = 50GeV (right). The dashed curves show the optimum configurations for the standard optimization 
(shaded region from the lower right panel of Fig. \E^. True parameter values of siia^ 29^^^ = 10~^ and 
5qp^ = 37r/2 have been assumed. The diamonds show the optimal configurations for non-standard physics, 
whereas the circles correspond to the IDS-NF standard choices Li = 4 000 km and L2 = 7 500 km. 
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7 Summary and discussion 

We have discussed the optimization of a neutrino factory for non-standard interactions (NSI) 
in the neutrino propagation in terms of muon energy, basehnes, and oscillation channels. 
Our study has been based on both analytical formulas, and a full simulation of the IDS-NF 
(international design study of a neutrino factory) baseline setup with GLoBES. We have 
considered all possible non-standard parameters e^n, and have also included the complex 
phases of the off-diagonal elements. 

As far as the different e™^ and different oscillation channels are concerned, we have identified 
the v^ appearance channel as the main contribution to the |e™ | sensitivity and the i/^ disap- 
pearance channel as the main contribution to the |e"^| and |e!^| sensitivities. Furthermore, 
|e^| and |e™ | cannot significantly be constrained beyond the current bounds at the neutrino 
factory (c/.. Table [1]), and the |e^| sensitivity can be directly related to the matter den- 
sity precision measurement dominated by the v^ appearance channel. Therefore, we have 
focused on e^, e™ , and e™ in the main line of this study. For example, we have presented 
analytical formulas for these quantities and the corresponding channels. Note that we have 
only considered one non-standard parameter at a time, because we have demonstrated that 
the two-parameter correlations are either unimportant if appearance and disappearance in- 
formation is used, such as in the e™-e^ or ej^^-e^ planes (or e^-e^^, e^-e^ planes), they can 
be analytically understood in a straightforward way, such as in the ej^^-e™ plane, or they 
can be related to the matter density uncertainty, such as in the e^-e^ plane. 

We have also considered the silver v^. -^ v^ channel for non-standard interactions, and we 
have only found a synergistic, but small contribution to the |e^| sensitivity if E^ > 25 GeV. 
This finding is in tension with the current IDS-NF baseline setup: At least for NSI, if the 
muon energy is chosen to be as low as 25 GeV, the silver channel will be hardly useful. In 
combination with the standard oscillation parameter optimization from Ref. [9], we conclude 
that the tension can only be released if either the muon energy is increased, or if the emulsion 
cloud chamber is removed from the IDS-NF setup. Except for the silver channel contribution, 
we have demonstrated that the NSI sensitivities do not significantly improve anymore as 
a function of the muon energy if E^ > 25 GeV, unless e^ is searched for at intermediate 
sin^ 26'*3^*' ~ 0.001 in well-known matter density environments. 

Furthermore, we have revisited the optimization for the standard oscillation parameters 
as a function of the two baselines of the two main detectors, and we have found that the 
optimal detector locations are consistent with the IDS-NF setup. We have then established 
the robustness of this optimization with respect to a possible NSI pollution, even though 
the absolute sensitivities become deteriorated. As the next step, we have studied the NSI 
sensitivities as a function of the two baselines. We have found that the optimization of the 
standard oscillation parameters is consistent with the one for |e^|, while for |ej^| and |e^|, 
even longer baselines are in principle preferred (basically, as long as possible within the 
Earth's diameter). Note that in all cases, one very long baseline (> 7000 km) has turned 
out to be a key component for the non-standard matter effect measurements. 

Our main results are summarized in Fig. [71 Obviously, E^ = 25 GeV provides an excellent 
sensitivity to all standard and non-standard performance indicators, while lower energies 
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Figure 7: Summary for the optimization of a neutrino factory as a function of the muon energy. The 
dark bars represent E^ — 50 GeV, the medium hght bars E^ ~ 25 GeV, and the hght bars E^ — 5 GeV. The 
upper group of bars represents the standard optimization (in terms of the sin^ 20i3 reach), the middle group 
of bars represents the standard optimization (in terms of the sin^ 26*13 reach) including e™ marginalized 
over, and the lower group the non-standard optimization (in terms of the |e™^| sensitivity). Here the IDS-NF 
setup is used with two baselines at 4 000 km and 7 500 km. Both the sensitivities without silver channel, 
as well as with an advanced silver channel detector Silver* arc shown in all cases. As a benchmark point, 
sin^ 26*13"° — 0.001 and S^p'^ — 37r/2 has been chosen, as well as a true normal hierarchy. 



23 



Performance indicator 



90% C.L. 



3a C.L. 



5(7 C.L. 



Standard oscillation physics 

sin^ 2^13 

Normal hierarchy (for 6^^^ = 37r/2) 

Maximal CPV (5^^^ = 37r/2, NH) 



4.25- 10-5 
2.27- 10-5 
1.49- 10-5 



1.22 ■ 10-^ 
5.93 ■ 10-5 
4.68 • 10-5 



3.03-10-4 
1.23-10-4 
1.23-10-4 



Standard oscillation physics polluted by non-standard e^ 

^2 2^13 8.13-10-5 2.04-10-4 

4.00-10-5 1.01-10-4 

4.69-10-5 1.39-10-4 



sm ZC713 

Normal hierarchy (for 6qp'^ = 3tt/2) 

Maximal CPV ((5g^^ = 37r/2, NH) 



4.88 - 10-4 
2.29 - 10-4 
5.52-10-4 



Non-standard oscillation physics with real e 

^Tu (with Ime™ 



m 
al3 



"e/i 



" (with Im e™ : 
"^ (with Ime™ 



0) 

0) 

= 0) 



-fMIU, 



Non-standard oscillation physics with complex e^ 



[-1.77 


10-^ 


-3.70 


10-3, 


-6.33 


10-3, 


1.71 


10-3] 


3.26 


10-3] 


5.98 


10-3] 


[-4.46 


10-3, 


-9.18 


10-3, 


-1.37 


10-2, 


3.51 


10-3] 


5.98 


10-3] 


0.93 


10-2] 


[-3.69 


10-4, 


-6.76 


10-4, 


-1.16 


10-3, 


3.68 


10-4] 


6.74 


10-4] 


1.15 


10-3] 


[-1.37 


io-\ 


-2.77 


io-\ 


-3.48 


10-1, 


1.23 


10-1] 


2.26 


10-1] 


3.86 


10-1] 


[-1.90 


10-2, 


-2.64 


10-2, 


-3.58 


10-2, 


1.89 


10-2] 


2.59 


10-2] 


3.55 


10-2] 


[-1.90 


10-2, 


-2.62 


10-2, 


-3.57 


10-2, 


1.90 


10-2] 


2.62 


10-2] 


3.57 


10-2] 



~efi\ 



3.41 - 10-3 
4.74 - 10-3 
1.80- 10-2 



5.71 - 10-3 
9.36 - 10-3 
2.22 - 10-2 



8.08-10-3 
1.75-10-2 
3.31 - 10-2 



Table 1: Summary of the sensitivities achievable in the IDS baseline setup (Li = 4000 km, L2 = 7500 km, 
Efj^ = 25 GeV) without a silver channel detector for different CL (1 d.o.f.). The first two groups show the 
sin^ 26*13 reaches, the last two groups the e™^ reaches. In the case of real e™o, we give the positive and 
negative limits separately. Here sin^ 26*'™'' = 0.001 and S^p° ~ 37r/2, as well as a normal true hierarchy 
have been used as a benchmark point. Note that especially the e™ sensitivity depends on this benchmark 
point. 
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(represented by the light bars) are unfavorable for both standard and non-standard per- 
formance indicators. Higher muon energies or the addition of a silver channel detector at 
4000 km do not yield any significant improvement, neither for standard oscillations, nor for 
NSI. Furthermore, we show in Table [T] the expected sensitivities for the IDS-NF baseline 
setup 1.0, i.e., E^ = 25 GeV, Li = 4 000 km, L2 = 7500 km, but without the silver channel. 
For the standard oscillation parameters, the sensitivities become somewhat worse if there is 
an NSI pollution from e^ (second group versus first group of sensitivities), but the orders 
of magnitude do not change. In addition, we read off the table that for |e^| and |ej^|, the 
neutrino factory will hardly improve the current limits (see Ref. [2] and references therein). 
For |e^|, the current limit can be improved by about two orders of magnitude, for ej^, 
about one order of magnitude (if 0J^ is assumed to be free), and for |e™ |, about two orders 
of magnitude. For e^, the improvement depends very much on sin^ 26'i3 and S^p^, as well 
as the matter density uncertainty, as we have discussed in Sec. El In general, a factor of 
a few may be expected. We believe that this overall performance is very impressive, but 
remember that the current bounds might be improved by the time a neutrino factory is 
actually built. 

Of course, our study has been based on a particular type of new physics, namely non- 
standard neutral current interactions affecting the neutrino propagation. One may ask 
the question whether, within the IDS-NF baseline setup, our results can be qualitatively 
generalized to other types of new physics. In particular, are there scenarios which constitute 
a physics case for u^- detection? In general, there might be two qualitatively different 
approaches to search for new physics using the information from Ur events: 

1. Use the spectral dependence of the z/^ events and test its consistency with standard 
oscillations, and with new physics scenarios. This is what we have done in this study 
for the specific case of non-standard interactions in the u^. -^ v^- (silver) channel. 

2. Test unitarity by using information from all flavors. 

For approach 1, the v^- appearance channels can only contribute significantly in the unlikely 
case that the impact of the new physics is much larger in these channels than in the others. 
After all, the golden and disappearance channels provide much larger statistics, and will 
therefore typically dominate the measurement, if appropriate baselines and neutrino ener- 
gies are used. Note that the z/^ — >■ Vj. channel may not be feasible at all because the high 
event rates in this channel might prohibit successful reconstruction. For approach 2, one 
could either use neutral currents as a signal (which would not require a dedicated Vr de- 
tector), or consider the flavor sum of charged current event rates. For the neutral currents, 
systematical uncertainties and the charged current contamination will limit the measure- 
ment to a precision of, perhaps, a few percent [78]. For the charged currents at the neutrino 
factory, the weakest link will probably be the detection of electron neutrinos (preferably 
with charge identification), not the detection of v^-. Electron neutrino events are very dif- 
ficult to reconstruct using an iron calorimeter, because electrons produce electromagnetic 
showers. The associated uncertainty in the event rates is expected to be of the order of a 
few per cent, so that summing the charged current events over all flavors will not yield a 
sensitivity significantly exceeding that from neutral currents. Therefore, we expect that the 
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conclusions from this study concerning the silver channel are likely to be translated to many 
other new physics cases as well, which, however, needs to be proven in specific studies. 
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